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On quarkonium in an anisotropic quark gluon plasma 



We reconsider a plasma with an anisotropy imposed on the momentum distribution of the 
system and study the real time static potential for quarkonia. The distribution function is 
normalised so as to preserve the particle number in an ideal gas, as required in the Keldysh- 
Schwinger formalism. In contrast to recent findings without this normalisation, a weak 
anisotropy does not lead to an increase in the melting temperature for bound states. To test 
for the maximal effect, we also investigate a gluonic medium in the limit of an asymptotically 
strong anisotropy. The spectral function of quarkonium is calculated for this case and found 
to be in remarkable agreement with the corresponding results for an isotropic medium. 
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Abstract 



1 Introduction 



The QCD-plasma with an anisotropic momentum distribution has been subject to several 
recent investigations as a simple model for the early stages of heavy ion collisions [H El |3] . 
Such a system is obtained by replacing the isotropic momentum space distribution functions 

by 

fik,0 = N{0 /iso [Vk' + ^k-ny) , (1) 

which removes modes with momentum components along the direction of the anisotropy, n. 
Here, ^ parametrises the strength of the anisotropy and N{^) is a normalisation factor with 
-^(^ = 0) = 1. Interest has particularly focused on the question how such an anisotropy 
influences the dissociation of non-relativistic quarkonium systems within potential model 
studies in O [6]. The general approach in those works was to consider small anisotropics 
^ < 1 and to calculate the correction to the propagators and the static potential to first 
order in ^. The effect of the anisotropy on the normalisation was neglected, i.e. N{^) = 1 in 
those works, and with the exception of [7] all employ the equilibrium Kubo-Martin-Schwinger 
(KMS) condition, which strictly speaking does not hold in the anisotropic case. The general 
conclusion of those investigations was that the anisotropy tends to decrease the effect of 
Landau damping and thus to increase the dissociation temperature. 

In this paper we assess the extent to which these results are modified by taking a ^- 
dependent normalisation factor into account and solving again the Dyson- Schwinger equa- 
tions for the non-equilibrium system. This is motivated by the fact that modifying only the 
argument of the isotropic distribution function in Eq. ([T]) simply removes particles/modes 
with momentum components along the anisotropy direction, i.e. the plasma gets diluted with 
growing ^. The momentum distribution defines the thermal parts of the bare propagators 
in the real-time formalism and thus parametrises a gas of non-interacting particles. Hence, 
the anisotropy must not affect the average number density of particles. This results in a 
normalisation condition for all ^ [HE], 

» = /|^/..o(p) = /|^/(p.e). (2) 

with the solution iV(^) = a/1 -|- ^. For example, the integral over the Bose distribution func- 
tion, /iso(p) = nsip), evaluates to n = 3.606 T^/(27r)^. Using iV(^) = 1 for ^ 7^ gives 
instead n = {2.55, 1.09, 0.3} TV(27r)2 for ^ = {1,5,100}, respectively. Even for "small" 
anisotropy ^ = 1, the dilution amounts to a significant 30% and affects observables. We 
find that the modified normalisation largely compensates the effects observed previously, 
such that for ^ = 1 there is little difference in the static potential between the isotropic 
and anisotropic cases. We then proceed to study the 'maximal anisotropy effect' by consid- 
ering the limit ^ — > 00, in which the plasma is confined to a thermal plane in momentum 
space perpendicular to the symmetry axis of the system. Somewhat surprisingly, also in 
this case the quarkonium spectral function is qualitatively in agreement with the isotropic 
situation. For simplicity a purely gluonic system is considered, the inclusion of light fermions 
is straightforward and will merely introduce a modification to self energy prefactors |9|. 
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2 The static potential, propagator and self energy 



For our calculations we use the Keldysli-Schwinger real time formalism jTOj, similar to 
There one considers matrices of correlation functions, Cjj, with indices i,j G {1, 2} denoting 
the forward and backward portions of the Keldysh contour. The real time static potential 
in a thermal medium may be defined through the quarkonium Wightman function in the 
infinite quark mass limit, i.e. the Wilson loop in Minkowski time evaluated in a thermal bath 
pT] , whose time evolution is governed by a 'Schrodinger equation': 



N' 



-Tr < - 



w > , 



idtC2i{t,r) = V{t,r)C2i{t,r) 



(3) 



The static potential corresponds to the large time limit of this evolution, where the Wightman 
function may be replaced by Cu, 



Vir) = lim V(t,r) = lim — — — 

t^oo t^oo Cii{t,r) 



(4) 



and to leading order in perturbation theory is equivalent to a matching coefficient appearing 
in a non-relativistic effective theory for finite temperatures [121 UHl El]- Calculating the 
Feynman diagrams corresponding to the Wilson loop through order g'^ in the real time 
formalism, one finds the generic expression for the potential in an arbitrary (non-)equilibrium 
system [T5] . 



V{r) 



(2^ 
(27r)3 



1^. 



'1 



'ikr\ 



(5) 



The tilde indicates a HTL-resummation of the gluon propagator. Thus we need to evaluate 
the resummed retarded and symmetric propagators for the various degrees of anisotropy we 
wish to discuss. The propagators are solutions to the Dyson-Schwinger equations. 



Gb. — <JR + Gr ■ IIr ■ Gr, 

Gs = Gs + Gr ■ IIr ■ Gs + Gs • 11^ ■ Ga + Gr ■ lis ■ Ga 



(6) 



The bold face notation indicates matrices in Lorentz and colour space. In the following all 
colour matrices are ~ 6ab and will be suppressed to simplify the notation. The retarded 
self energy for an anisotropic plasma has been calculated in HTL-approximation in [16], and 
the retarded propagator in a general covariant gauge in [1]. The resummed retarded gluon 
propagator in a general covariant gauge can be written in the form 



G 



4 2 



(7) 
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Here the symmetry vector n'^ = (0, n) represents the anisotropy, K = {u, k), X is the gauge 
parameter and the electric and magnetic propagators are 

A^^ = (K' - a~j){u;' - (3) - 6'[k' - {n - Kf] , 

A^^ = K^-a . (8) 

Exphcit expressions for the tensors A — D and the structure functions a{K,^) — 6{K, ^) for 
arbitrary anisotropics are hsted in Eqs. (2-5) in |1] and Eqs. (B1-B4) in [T], respectively. 
Note however, that the structure functions given there are for N{C,) = 1 and need to be 
multiplied by N{^) = y/1 + ^ for our purposes. Since we are mostly interested in applications 
for potential models, it is expedient to take the static limit u; ^ already at this stage to 
facilitate the calculations. In this case the bare symmetric propagator vanishes and, using 
the fact that the retarded and advanced propagators are related via = {G^^)* in a 
general homogeneous system [10], the second of Eqs. ^ simplifies to (see also [9]) 

Gs = ■ Hs ■ (9) 

The propagator in the static limit is (with the notation k"^ = k'^, cos 6k = n ■ k/k) 



G>^'iu = 0,k) = AA[C"^"-An+AG 

A 



k^ 



l^k' + mi + mt^B^^" - {k' + m')C^'' + i—^D 
' sm t>: 

K'^K" (10) 



with now 

^l^ = e + ml, Ag' = (P + m^ +m2)(fc2 + m2) -m^ (11) 
The tensor basis is simplified and takes the form 



A'^^'{k) = A^^{k) = {] A'^{k) = 6ij-k'y/k^ B^^^ik) = S^oS^o 
C^i'{k) = C^ik) = C'^{k) = feh^ D^'^ik) = d^oh" + 5^on 



(12) 



with = A'HK A rescaling of the tensors ^B^"" B^"" , ^D^"" D^^, as well as 



,,2 



A;2 



Ag Ag, — Ayi — * Aa has been introduced with respect to Eqs. (7 8) to keep all 
components of the retarded propagator explicitly finite in the static limit. The effective 
masses are the static limits of the structure functions [1], 



m'i = lim a , m'i = lim -3 , ml = lim 7 , = lim lm(S) . (13) 



3 Weak anisotropy, ^ < 1 

In order to discuss the effect of a weak anisotropy, we employ the distribution function 
Eq. ([T|, including the non-trivial normalisation with iV(^), and expand the resummed prop- 
agators up to linear order in ^. For the longitudinal component of the retarded propagator 
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appearing in Eq. ([s]) we then have 



-|- 777,2 _|_ ^2 ^^^2 _|_ ^2^ _ ^4 



(14) 



All ^-dependence resides in the effective masses. To linear order the dimensionless combina- 
tions ml = ml/m'jj read 



ml = cos^ 9k, ml = I + ^{cos^ 9k - ]:] 
3 ^ o 



m^ 



sin 9k 



IT 



mx 



^ — sin 9k cos 9k ■ 



(15) 
(16) 



Note that differs from the expression given in [T] due to the normalisation factor N{^) = 
v^l + ^ in the distribution function Eq. (LLI). The retarded propagator to linear order then is 



r<oo 



1 



m 



D 



'cos^ 9k - 



6' 



(17) 



Expanding also the static potential, V{r) = Vo(r) + ^Vi{r) + . . ., we find for the correction 
of the real part due to the anisotropy, 



(1 - cos(fc • r)) 



1 



{k^ + miy 



cos^ 9i 



For the imaginary part we need in addition the symmetric propagator and the symmetric 
self-energy, cf. Eq. ([9]). In the static limit and for soft external momenta fcj -C Pi the latter 
takes the form [9] 



= Svr/ivl J + fiP + kMv, ■ v^) . 

Straightforward calculation gives the following expansion for the symmetric propagator 



(19) 



2n 



m 



— 1 



D 



k{k^ + miy (3 
1 



1 + ^ 



7r^-3C(3) 3.2/ 

sm 9k 

4 



4 



/5 k{k^ + miy 



fcos^ 9i, 



1, 



(20) 



where is the Riemann zeta function. For N{C,) = 1 the first fraction in square brackets is 
absent and 1/6 — >■ 2/3 in the last line. In this case we reproduce the symmetric propagator 
from [7] . The leading correction to the imaginary part of the potential due to the anisotropy 
is given by 

— il-,osk.r)--Gfl^ (21) 
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Figure 1: Real (left) and imaginary (right) part of the static potential. The blue curve 
corresponds to the isotropic plasma, = 0, the other curves to ^ = 1 with normalisation 
factor N{^) = 1 (orange) and N{^) = y/1 + ^ (red). The potential is evaluated for 9r = ^ 
(solid) and 9r = (dashed). 

The remaining integrations can be done numerically, and the results for the real and imag- 
inary parts of the static potential are shown in Fig. [T| where we have defined the angle 
between r and the anisotropy axis, cos 6r = r ■ n/r. For the normalisation N{C,) = 1 it 
has been observed that an increasing anisotropy lowers the real part of the static potential 
towards a Coulomb potential ^ . However, this follows trivially for ^ — cx) since the medium 
is disappearing in this case, and hence tjid ^ 0. By contrast, the effect of the anisotropy is 
much weaker when the normalisation iV(^) = a/1 + ^ is chosen. Note also, that in this case 
the potential is even slightly enhanced at small distances. The difference between the two 
normalisations is even more pronounced in the imaginary part, which is caused by Landau 
damping and introduces a broadening of the spectral function [17] . A weaker damping effect 
is seen for N{C,) = 1, leading to the conclusion that bound states should melt at higher tem- 
peratures [51IS1E]- For N{^) = y/C, + 1, this effect is modified and persists only at very large 
distances. For smaller distances, on the scale of the bound states, there is only very little 
difference to the isotropic case. Uncertainties due to non-perturbative corrections to the dif- 
fusive physics underlying the existence of the imaginary part were found to be significantly 
larger in recent studies [T8| [19]. 

4 Strong anisotropy, ^ ^ 00 

Having seen that, with a distribution function normalised to maintain particle number in an 
ideal gas, the effects of a weak momentum space anisotropy on the static potential are rather 
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mr)/k m£)/k 

Figure 2: The effective masses of stable electric and magnetic modes (shown on the left for 
^fc = 7r/4 and k < mi)) indicate a screening of long range interactions. For 6k < vr/4 the 
gluon propagator has poles on the imaginary u-axis due to the Weibel instability. The growth 
rate 7 = Im(u;) of the associated unstable modes (right) increases towards the symmetry 
axis. 



small, we are now asking how a maximally anisotropic plasma would behave. We begin by 
returning to the gluon propagator in a general covariant gauge, Eq. ([T]), and discuss some 
aspects of the dispersion relations which have been derived in [2]. 

4.1 Dispersion relations 

In the limit of maximal anisotropy, the structure functions a{K) — S{K) take on an analytic 
form [2], facilitating a numerical evaluation of the dispersion relations of the retarded gluon 
propagator, = and = 0. The associated modes are referred to as electric and 
magnetic modes respectively (for an analysis of dispersion relations at arbitrary anisotropics 
see [U E]). They can be conveniently expressed by introducing effective pole masses m(fc), 
which satisfy the relation 

cj2 = p + m2(fc) (22) 

at poles and branch cuts of the electric and magnetic propagators. The results from a 
numerical determination of the pole masses are illustrated in Fig. [2] and summarised in the 
following. The dynamical screening masses of the lowest electric mode as well as the magnetic 
mode are found to be comparable in magnitude to the screening masses in an isotropic 
medium with a moderate angular dependence of both masses towards higher momenta. 
The findings indicate that strongly coupled interactions over large distances are suppressed 
leading to a deconfinement of colour charges. This is a prerequisite for the possibility of a 
consistent perturbative treatment using the HTL approach. 
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For 9]^ < tt/4 the gluon propagator has additional poles along the imaginary u-axis indi- 
cating the existence of unstable field modes. The instability, which is present in a general 
anisotropic QCD plasma, is speculated to be responsible for the rapid isotropisation of the 
quark gluon plasma observed in heavy ion collisions [20]. The growth rate 7 = Im(ci;) of un- 
stable modes is observed to increase and to peak at higher momenta towards the symmetry 
axis. 



4.2 The static potential 

Next, we take the static limit of the maximally anisotropic propagator in order to discuss 
the static potential. At maximal anisotropy the appropriate effective masses read: 



m ^ = 

4tan^^ 



nmj-) 2 7rm|, 2 7rm|, / ^ \ 2 7rm|, cos6'fc 

^P = ^^ '^^ = — [^ + ^^J^ '^^= 4 sin^^, 



k 

(23) 

We also need to re-evaluate the symmetric self energy, Eq. (19). Since the momentum 
distribution function for the anisotropic medium /(p) becomes sharply peaked around Vp ■ 



n = upon sending ^ to infinity, the following identity may be derived 

/oo 
dx n-BipVl + x^) = 6{vp ■ n)h{p). (24) 

The self energy thus takes the following analytically tractable form: 

inr(fc) = / « Kp)Kvr> ■ n)6{Vp ■ v,) [l + Mp)^±^5(^;, ■ n)) . (25) 

The Bose enhanced part of the self energy is amplified by a factor \p + k\/k p/k which will 
translate into an amplification by T/fc, indicating that scattering processes are dominated 
by Bose enhanced scattering into the thermal plane. The self energy is diagonal, 

= ^ = w ^ ■ ">) • = = " ■ '''' 

with m|) = y^f^T^ and the prefactors Cx are obtained from the following integrals 



/ dpp^h{p) = 2.9231 



TT 



fi /-oo 

c± = —(3^ dpp^h^p) =3A230. 

TT^ Jo 



(27) 



The symmetric propagator can now be calculated from (|9]) in a straightforward fashion. Its 
longitudinal component can be separated in contributions dominated by external momenta 
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Figure 3: The real part of the anisotropic static potential, normalized to V{oo) = and the 
imaginary part shown as a function of the radial separation r for = 3 and = 3. The 
cylindrical symmetry of the anisotropic imaginary part is illustrated by varying Or- 



perpendicular and parallel to the symmetry axis 
iG°°(fc) = iG^(fc)+i4(fc) 



c±ml) 



lUg[K) - , ™2 , ™2V;.2 , ™2^ ^412 ' ^"^^ 



The first component corresponds to Bose enhanced scattering into the thermal plane and 
is enhanced in |^ ^ 1 while the second component, corresponding to scattering away from 
the thermal plane, is suppressed even for sin ^ due to the diverging effective mass 
in this limit. Contributions from the symmetric propagator will therefore be dominated by 
momenta perpendicular to the symmetry axis, 

iG°°(fc)^iG^(fc). (29) 

We are now ready to discuss the static potential in a maximally anisotropic medium, by 
inserting these ingredients into Eq. ([s]). Upon choosing a spherical coordinate system with 
the anisotropy vector n pointing in the z-direction, cos 6r = r-n/r and with the abbreviation 
X = cosdk, one angular integration can be performed immediately leading to 

^1 — Jo{krV 1 — x'^ sin 6^.) cos {krx cos 6*^) j 



k'^{k'^ + + m. 



(/c^ + + m'^){k'^ + m^) — mj 



^' (30) 



4 ' 
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where Jn is the Bessel function of the first kind. The expression is readily evaluated numer- 
ically using standard integration routines p2]- The potential contains a linear divergence 
which we subtract by normalising V{\r\ —* oo) = 0. The resulting real part displays an 
approximate spherical symmetry. The radial profile, as depicted in Fig. [3] for = 3, turns 
out to be in good agreement with the standard Debye screened potential. For the imaginary 
part we obtain in a similar way 

^-'nr) = /dMl - Mkr.inO.))^^,^L-^^ . (31) 

The imaginary part thus only depends on the separation rsinOr of the static quark pair 
in the thermal plane perpendicular to the anisotropy vector. Unlike the previously studied 
quantities the imaginary part turns out to be quantitatively different from its isotropic coun- 
terpart [11]. The damping of the static quark correlator is significantly weakened towards 
the anisotropy axis for reasonable values of the coupling. 

4.3 Quarkonium spectral function 

In order to compare a more physical quantity between the isotropic and anisotropic cases, 
we also calculate the quarkonium spectral function. Since this has not been done before for 
anisotropics, we present it in some detail. For a qq system with finite constituent mass M 
the following Schrodinger equation must be solved |llj . 



idtC2iit,r) = (^2M-^ + Vir)^ C^2i(t,r), iC2i(0,r) = -6N5^^\r) (32) 

where V{r) is the real time static potential discussed in the previous sections. For r = the 
quarkonium correlator coincides with the quark current correlator C2i(t, 0) = ( J'^(t) J^(0)), 
where = ip'j^ip- The quark pair is assumed to be created at time t = and the initial 



condition for the Schrodinger equation (32) ensues from the afore mentioned relation. The 
spectral function for quarkonium at rest p{uj) can be calculated from the difference of the 
two off-diagonal elements of the heavy-quarkonium correlator, 

Pico) = Jdt e--* {iC2i(t, 0) - iCi2(t, 0)} . (33) 

In thermal equilibrium both components are related via the Kubo- Martin- Schwinger con- 
dition, which does not hold for the system under consideration. In the operator formalism 
the difference between the two correlators takes the form 

C2i{t,0)-Cu{t,0) = ^Tr {<7j^(t,0)>(0,0)} - ^Tr |<tJ^(0, 0)>(t, 0)} , (34) 

where the trace is over all energy eigenstates and a is the statistical operator of the anisotropic 
system. Using the cyclicity of the trace, C12 is expressed as 

Cu{t,0) = |Tr{j^(t,0)a>(0,0)} . (35) 
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The difference between both contributions is that in C12 the density matrix acts on the 
Hilbert space of the plasma with the quark pair present. For 2(3M ^ 1, where the quarko- 
nium threshold ~ 2M is much higher than the average particle momentum f3~^ in the thermal 
plane, the contribution of the correlator C12 will be substantially suppressed. This justifies 
the following approximation which is also employed in the equilibrium case [17| . 



p{uj) ^ J dte-"^'C2i{t,0). (36) 

The Wightman propagator C2i{t, r) for a static quark pair is obtained at large positive times 
as a solution of the Schrodinger Eq. (|4|. Similarly, using iCij{—t, r) = [iCij]*{t, r), the static 
potential in the infinite past is 

V\r) = hm '^f'f'':^ = -/Cf / ^(1 - cosfc ■ r)G°° = V*ir) . (37) 



The spectral function can thus be computed from the relation 

POO 

p{uj) = 2 dt {cos (c<;t)Re(iC2i) - sin {ujt)lm{iC2i)} . (38) 
Jo 

The numerical solution of the Schrodinger equation for the Wightman correlator C21 is non- 
trivial. The task consists in the interesting problem to solve a parabolic complex differential 
equation in three dimensions with a potential consisting of a parametrically large spherically 
symmetric real part and a small cylindrically symmetric imaginary part. Due to the positive 
parity and the symmetry of the potential with respect to the anisotropy axis, the correlator 
may be expanded in a subset of spherical harmonics, 

C2i{t,r) = Q(t,r)Y°(^^r), 

V{r) = \/(-)(r) + iV^(^")(r)YO(^^.), (39) 

with V^^'^^r) and Vi^^"^\r) denoting the real and imaginary parts of the potential respectively. 
With the rescaling 

Q(^,r) = ^^e-2*^ (40) 

the one-dimensional Schrodinger equations for the components of the spherical expansion 
take a numerically tractable form, 

= -j^{^2- ^^^) + y'--'-^ + ^-PlmnV^-'-n • (41) 

The tensor Vimn is defined by a unit sphere integration over a product of three spherical 
harmonics which can be evaluated exactly via Wigner-3j symbols [23] : 
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Without the couphng between different spherical components, a relation of the following 
type needs to be evaluated for every timestep dt, 

ui{t + dt)=e-'^'^%{t), (43) 

where the differential operator H/ is given by 

To ensure a stable time evolution we employ the Crank-Nicolson scheme by introducing 
the second order approximation 

1 + ^iHdt^ ui{t + dt) = - ^iHdt^ ui{t) , (45) 

which is solved as a tridiagonal matrix equation for the discretized system. The full problem 
is solved without sacrificing the stability of the algorithm by treating the parametrically 
small coupling between the various spherical components as a linear perturbation, 

ui{t + dt) = ui{t + dt) + dt 5^ VimnVj;!"'^ n„(t) . (46) 

The system is discretised on a finite lattice with spacing a and lattice sites [1, . . . , N^], 
choosing the boundary conditions 

ui{0) = and ui{Nl + 1) = (47) 

for all spherical components with / limited to the range / G [0, . . . , Imax]- The initial condition 



(32) is discretised as follows [T7] . 



with r = na being the discretised radial coordinate. All other spherical components m;^o(0, r) 
are set to zero. The scale for the numerical determination of the bottomonium spectral 
function for = 3 shown in Fig. [4] is set by the QCD scale Aj^g ~ SOOMeV in the 



MS scheme. The bottom quark mass M = 14.15A]g;5r is chosen in agreement with [T7] . 
Introducing a spherical cutoff Imax = 10 the system of coupled one dimensional Schrodinger 
equations is solved on three different grids G {2000,4000,8000} with a fixed physical 
size of aN^ = 400/A]g;g^. A time-step of dt = 0.02a is chosen and in all three cases the 
system is evolved for a fixed time t = 265/AjYg. The spectral function depends linearly on 



the lattice spacing a since the coupling between the different spherical components of (41 ) is 
introduced as a linear perturbation. To take into account a^-corrections, the spectral function 
is extrapolated to the continuum by fitting the a-dependence of the spectral function at a 
given frequency against a quadratic polynomial. The melting of the quarkonium resonance 
due to Landau damping effects is in qualitative agreement with the decay of the resonance 
in thermal equilibrium obtained using the same approach [T7] . 
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Figure 4: Left: The melting of the bottomonium spectral function in a maximally anisotropic 
quark gluon plasma at a fixed coupling of = 3 (a^ ~ 0.25). Right: The polynomial 
extrapolation of the spectral function using three different grid spacings. 

5 Conclusions 

We have reconsidered the effects of an anisotropic momentum distribution on the real-time 
static potential of a heavy quark pair in a plasma using a normalised distribution function 
which keeps particle numbers constant in the ideal gas limit. At small anisotropics, previ- 
ously observed effects were shown to be mostly due to a particle dilution stipulated by a 
different normalisation factor. In order to probe for a larger effect we then looked at max- 
imal anisotropy. The calculation of the quarkonium spectral function indicates a melting 
of the quarkonium resonance similar to the one for an isotropic medium due to a screening 
of the Coulomb potential and a dissociation of bound states by Landau damping effects in 
the thermal plane. Strongly coupled interactions over large distances are suppressed due to 
the generation of nearly isotropic screening masses comparable in magnitude to their equi- 
librium counterparts. For small angles with the anisotropy axis, the gluon propagator has 
additional poles along the imaginary cj-axis indicating the existence of unstable field modes. 
The growth rate 7 = Im{uj) of unstable modes is observed to increase and to peak at higher 
momenta towards the symmetry axis. 

The overall physical picture that emerges from a weak-coupling analysis of this system 
indicates that the strongly anisotropic state of matter created in the initial stages of heavy-ion 
collisions is characterised by an early deconfinement of colour charges and a rapid diffusive 
transport of these charges in the collision plane. The most striking result is that observables 
related to quarkonium are very similar to their counterparts in an isotropic medium, making 
it difficult to distinguish a highly anisotropic plasma from its isotropic counterpart on the 
basis of heavy quark observables. 
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